      function esat (t)
c        tensao de saturacao
      implicit double precision (a-h,o-z)
      esat=618.78*exp(17.269*(t-273.16)/(t-35.86))
      end